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Abstract 

We study the cosmological impact of tiie super symmetric DFSZ axion model. Extending recent 
works, we first provide a comprehensive analysis of thermal production of the DFSZ axino con- 
sidering all the possible scattering, decay and inverse decay processes depending on various mass 
parameters and the reheat temperature. Although it is hard for the DFSZ axino to be in thermal 
equilibrium, its coupling is still large enough to generate huge axino population which can turn into 
overabundant neutralino density. We examine the neutralino parameter space to identify the dark 
matter property depending on the Peccei-Quinn scale. As the Peccei-Quinn scale becomes higher 
resulting in longer axino lifetime, the neutralino dark matter appears in a lighter Higgsino-like 
LSP region or a more restricted Bino-like LSP region allowing a resonant annihilation through a 
CP-odd Higgs boson to meet stronger reannihilation. 
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I. INTRODUCTION 



The most attractive solution to the strong CP problem would be to introduce the axion, 
a pseudo-Goldstone boson of the Peccei-Quinn (PQ) symmetry which is supposed to be 
broken at vpq = 10^ — 10^^ GeV l|. The axion solution can be realized typically in two 
ways: (i) the KSVZ model assuming a heavy quark at an intermediate scale ^ fpg j2i; 
(ii) the DFSZ model extending the Higgs sector with PQ charged singlet fields j3|. In the 
supersymmetric extension of the axion models, the presence of the axino a, the fermionic 
partner of the axion, can change dramatically the property of dark matter. 

First of all, the axino can be the lightsest supersymmetric particle (LSP) in i?-parity 
conserving models. The axino in the KSVZ model decouples at a very high temperature 
and thus it must be very light {rria ~ 0.1 keV) to become a viable dark matter candidate j4|. 
In the gravity mediated supersymmetry breaking scheme, such a light axino mass could be 
arranged although it is generic to have the axino mass of order of the gravitino mass, 10^ — 10^ 
GeV |5|-[7|. In the DFSZ model or in the KSVZ model with lower reheat temperature, the 
axino may never be in thermal equilibrium, but thermal regeneration (the axino production 
from the scattering, decay and inverse decay of thermal particles) becomes an important 



source of the cosmic axino abundance 



14j . Although suppressed by the intermediate PQ 



symmetry breaking scale, vpQ, the axino interactions turn out to be efficient enough to 
generate a huge number of axinos unless the reheat temperature is much lower than the 
weak scale. 

If the axino has a typical mass in the range of 10^ — 10^ GeV and the reheat temperature is 
above the electroweak scale, the axino cannot be dark matter as its thermal population easily 
overdoses the Universe. Then, it has to decay into the LSP dark matter such as the usual 
neutralino LSP 



12 
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16| or gravitino [l7|. Since the axino lifetime is suppressed by Vpq, 
the axino decay may occur very late to produce the neutralino dark matter after its freeze-out 
causing overabundant neutralino density. This problem can be evaded if vpq takes its lower 
value to render the axino decay before the neutralino freeze-out, otherwise the neutralino 



LSP must have a large annihilation cross-section to allow a sufficient reannihilation 



15|- 



The thermal axino production in the KSVZ model has been studied extensively taking 
a simple assumption of the effective QCD interaction below the PQ scale Ulnl and more 



carefully considering a general situation 



13|. Axino couplings in the DFSZ model are dif- 
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ferent from those in the KSVZ model and their impact on the thermal axino production has 



been considered recently in Refs. 
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2|-|l4j. In this paper, we provide more general analysis 
of the DFSZ axino abundance including the scattering, decay and inverse decay processes. 
Then, we explore the Bino-Higgsino dark matter parameter region depending on the PQ 
symmetry breaking scale. 

The paper is organized as follows. In Section 2, we describe the supersymmetric DFSZ 
model to derive the dominant axino couplings. In Section 3, the DFSZ axino abundances 
coming from various channels are calculated and compared with each other. Since the 
DFSZ axinos are overproduced as far as the reheat temperature is above the weak scale, we 
investigate its impact on the neutralino dark matter property by taking some benchmark 
parameter points in Section 4. We conclude in Section 5. 



II. AXINO COUPLINGS IN THE SUPERSYMMTRIC DFSZ MODEL 

One of the hierarchy problems in the minimal supersymmetric standard model is the 
fi problem: "why the Higgs bilinear parameter, /i, sits at the TeV scale?". This problem 



can be nice 
Kim-Nilles 



y resolved in connection with the axion solution of the strong CP problem a la 
. Extending the Higgs sector with 



p2 

Wdfsz = ^HY^^uHd, (1) 

IVlp 

where Mp is the reduced Planck mass, one can introduce the PQ symmetry under which 
the Higgs doublet superfields Hu^d and the gauge singlet field P are charged. In order to 
break the PQ symmetry, one can introduce a PQ sector like 

WpQ = Xss(^PQ-'^^ (2) 

which will lead to (P) ~ (Q) ~ Vpq/V2. One can then get an appropriate /i term: = 
X_^{P)'^/Mp. This is nothing but the supersymmetric realization of the DFSZ axion model 
In this scheme, the axion supermultiplet is a combination of mostly P and Q, and a small 
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a. 



mixture of Hu and Hd [7[ • Thus, the axino interaction is given by the axino-Higgs-Higgsino: 

Cdfsz = CH-^a [HuHd + H^Hd] + h.c. (3) 

VPQ 

And after the electroweak symmetry breaking, the axino-top-stop couphng induced by the 
axino-Higgsino mixing also becomes important: 

C-DFSz,top = ct a[tt'' + W] + h.c. (4) 

VPQ 

where ch and q are order-one parameters depending on the PQ symmetry breaking sector 
and other couphngs proportional to smaller quark or lepton masses are neglected. As the 
axino couplings are of order /i/fpg < 10~®, the axinos cannot be in thermal equilibrium 
but can be produced thermally by the scattering, decay and inverse decay coming from the 
couplings ([3]) and (j4]). For our numerical calculation of the thermal axino production in the 
next section, we will put ch = 2 and q = 2, which corresponds to assigning the PQ charge 
1 to each Higgs supermultiplet, respectively, with the PQ sector Eq. 1^} Note that the 
axino-top-stop coupling arises only after the electroweak symmetry breaking and thus does 
not contribute to the thermal axino production if the temperature of the Universe is larger 
than the weak scale. 

Let us now compare the DFSZ axino couplings to the KSVZ axino couplings. In the 
supersymmetic KSVZ model, one introduces a heavy quark pair {X, X^) carrying the PQ 
charge: 

Wksvz = XxPXX' (5) 

from which the heavy quark gets an intermediate mass, Mx ~ Xx^pq/ V^- Below the scale 
Mx, the axino gets an effective QCD interaction: 



from which the previous calculation of the KSVZ axino abundance was preformed 



m 



^ Note that the physically relevant parameter is cn.t/vpQ so that CH,t also depends on the normalization 
convention on vpg. ct can be obtained from the axino-Higgsino mixing Eq. (|27|) with fi ^ mg. 
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FIG. 1: IPI axino-gaugino-gauge boson interaction at one-loop level for the DFSZ axino 

to show that the axino relic density is proportional to the reheat temperature Tr of the 
Universe. However, a recent study showed that such a T^-dependence does not occur if 
Tr is between vpq and Mx as the axino abundance is determined by the Yukawa coupling 



Xx ~ \/2Mx /vpQ and this property is app 
Let us recapitulate the result of Ref. 



ied to the DFSZ axino [13|.2 
13| to show that the DFSZ axino abundance is in- 
dependent of the reheat temperature. After the PQ symmetry breaking, we get the effective 
axino coupling with the Higgs and Higgsino fields ([3]) (which corresponds to the non-linear 
realization of the PQ symmetry with the parameter c = in Ref. 13|). This leads to the 
IPI axino-gaugino-gauge boson interaction of dimension 5 from the diagram in Fig. [H The 
corresponding amplitude is given by 



A 



9l 



_ -Ch I dx 

IQ'k'^\/2vpq Jo jo 



l-x 



dy 



- [pIx{1 -x)+ ply{l -y) + 2{pi ■ P2)xy] 



(7) 



for the axino and the SU(2) gauge supermulitiplet interaction. It is then easy to see the 
limiting behavior of the loop integral: if \p\ ^ fi, and if /i ^ \p\ assuming fi ~ mn, 
which tells us that the axino production through the effective axion-gaugino-gauge boson 
interaction is suppressed by /i^/T^ for T > /i if there is no other heavier gauge charged 
and PQ charged supermultiplet than the scale /i, which is the case for the "pure" DFSZ 
models.^ In this case, the axino production will be determined by the tree level interaction 
Eq. 03]) so that the DFSZ axino abundance becomes independent of Tp. Still if we assume 
some additional heavier gauge charged and PQ charged supermultiplet of some mass M 
larger than the scale fi, the axino production will be dominated by the effective interaction 



^ Note that a different claim was made in Ref. [3l ■ 

^ We refer the "pure" DFSZ to the models where gauge charged and PQ charged particles are only among 
the MSSM particles, and the "mixed" DFSZ to the models containing exotic gauge charged and PQ 
charged supermultiplcts in addition to the MSSM particles. 
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Eq. ([6]) between the energy scales Stt^/i and M as can be seen in the table 1 of 13|, and the 
axino abundance becomes proportional to Tr in this regime as in the KSVZ case. Previous 
studies on the DFSZ axino production with the Eq. ^ as the dominant interaction at high 
reheat temperature effectively correspond to such "mixed" DFSZ models. In this paper, we 
consider the "pure" DFSZ case with the Higgs supermultiplet as the heaviest gauge charged 
and PQ charged supermultiplet. So while the DFSZ axino abundance is independent of Tr, 
the KSVZ axino abundance is proportional to Tr and can be larger than the DFSZ axino 



abundance for larger Tb |12|. 



Note that the axino effective interactions with SU{2) x f/(l) gauge supermultiplets are 
obtained below the scale /i by integrating out the Higgs supermultiplets. In a similar way, the 
effective axino-gluino-gluon interaction ^ in the "pure" DFSZ model arises after integrating 
out the quark supermultiplet having the Yukawa coupling as in Eq. (jlj), along with its 
supersymmetric counterpart: the usual effective axion-gluon-gluon interaction solving the 
strong CP problem. 

Another distinction of the DFSZ axino from the KSVZ axino comes from the lifetime 
calculated from the couplings ([3]) and (Ej), respectively. As was shown in Ref. the 
thermal abundance of the KSVZ axino can dominate the energy density of the Universe 
for larger vpq. But the DFSZ axino lives shorter than the KSVZ axino and thus it always 
decays before it overdominates the Universe as will be discussed in Section [IVl 

III. THERMAL PRODUCTION OF THE DFSZ AXINO 

■"SZ axino has been computed previously considering the 



Thermal production of the D 
decay and inverse decay in Ref. 



121 ] and the scattering involving the axino-Higgsino-Higgs 



131 ] as depicted in Fig. [21 There 



coupling and gauge coupling of Higgs multiplet in Ref. 
are additional processes involving the top Yukawa coupling as depicted in Fig. [3] which can 
be potentially more important since the top Yukawa coupling is of order of unity. In this 
section, we combine all the possible processes to evaluate the axino yield for each channel. 
For the calculation, we take some benchmark points for the relevant particle masses including 
supersymmetry breaking effect both in the unbroken phase and in the broken phase of the 
electroweak symmetry. 

In Fig.m the thermal axino number density = n^/ s per the entropy density is analyzed 



6 



H 




a H 



H 



H W 




\ H 



(a) 



(b) 



FIG. 2: Examples of scattering diagrams (considered in Fig. 3-7 of Ref. [13|]) involving the axino- 
Higgs-Higgsino coupling and gauge couplings of the Higgs multiplet. There are also two other 
diagrams with supersymmetric gauge vertices and t, u-channel diagrams. All processes are included 
in the calculation. 
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FIG. 3: Examples of scattering diagrams involving the axino-Higgs-Higgsino coupling and top 
Yukawa coupling. There are also t, u-channel diagrams. All processes are included in the calcula- 
tion. 



for the benchmark points: fx = 1 TeV, = 100 GeV, 171^^ = 500 GeV,^ = 1-5 TeV, and 
rria = 100 GeV with vpq = 10^^ GeV. In this parameter space, the Higgsino decay H — )■ Ha is 
allowed to contribute to the thermal axino production in addition to the scattering processes 
from Fig. [2] and [3] (red line). As can be seen in the plot, three contributions are almost 
the same in the region of Tr > /i. In the slightly lower temperature, i.e. Tr < /i, the 
decay contribution is about an order of magnitude larger than the others because scattering 
processes are doubly suppressed by the Boltzmann factor while decay process has only one 
Boltzmann suppression. 

The second example is depicted in Fig. [S] for which we take the benchmark points: /i = 1 
TeV, niH^ = 100 GeV, mn^ = 500 GeV, m^- = 300 GeV, ma = 100 GeV with vpq = 
GeV. The only difference from the previous example is the stop mass. In this example, a 



and denote {rn^ + and (mf + where TO2 and mi arc soft mass terms of Hu and 



Hd, respectively 
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FIG. 4: The axino yield from the scattering diagrams in Fig. [2] (blue line), the scattering diagrams 
in Fig. [3] (red line) , and the decay H — >• Hd (green line) . 
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FIG. 5: Same as in Fig. Ubut with the parameter choice allowing a resonant axino production in 
the channel of Fig. El^b) 
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FIG. 6: Same as in Fig. but with parameter choice allowing the inverse decay HH 
line) 



a (purple 



small stop mass is considered to allow a resonant axino production in the channel t + 1 — )■ 
if — 7- a + if (Fig. [3l^b)). The resonant process can enhance the production cross section 
by factor of l/F^, where F is the decay rate of the intermediate particle, the Higgsino in 
this case. As the axino thermal production is obtained by integrating the cross section 
over T = [0,Tr], such a resonance effect is averaged out by the width of the resonance, F, 
resulting in the enhancement of the axino yield by factor of 1/F. In this example, we have 
a quite large decay rate of Higgsino, F ~ 100 GeV, leading to fi/T ~ 0(10) enhacement for 
the resonant production channel. Thus, the corresponding scattering process dominates the 
others by about factor 4. 

The third case, shown in Fig. [6l^ is for a heavy axino which can decay to a Higgs boson 
and a Higgsino: rria = 1 TeV and fi = 300 GeV with the other parameters kept same as in 
the previous case. As one can see, the axino yield from the inverse decay is more or less the 
same as the decay contribution in the previous examples. 

From the above calculation, one finds the axino yield Ya ~ 10^^(10^^ GeV/vpgY for each 



^ Our results in Figs. HHS] are different from Fig. 3 of Ref. lj| showing dependence as remarked in 
Section II. 
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FIG. 7: The axino yield at low reheat temperature below the electroweak scale with a decoupled 
benchmark spectrum 

allowed channel of scattering, decay or inverse decay as far as the reheat temperature above 
the electorweak scale (Tr > /i). In some cases, a resonant production can occur to enhance 
the axino yield in a certain scattering process. As a consequence, thermally produced axinos 
would severely overdose the Universe if they take a weak scale mass and are stable. 

The thermal axino production can be effectively suppressed if the reheat temperature 
of the Universe is below the weak scale. In such a case, the relevant processes for axino 
production strongly depend on a specific mass spectrum with broken electroweak symmetry.^ 
For the purpose of our presentation, we will discuss only a representative behavior taking a 
simplified benchmark point. To this end, we assume that the stop, Wino and CP-odd Higgs 
are heavy enough (~ 1.5 TeV) to be decoupled at low reheat temperature. Also we set 
rria = 500 GeV which is taken to be the reference value in the next section discussing heavy 
axino decay. For the Higgsino mass, we choose /i = 1 TeV allowing the decay H — )■ ha. 
Note also that the t ^ ta decay processes are possible. Fig. [7] shows the result for the 
axino yield in this benchmark point at low reheat temperature below 200 GeV. As we have 
observed in Fig. HI the decay contribution from the Higgsino is dominant over the scattering 

^ In this calculation, we use the zero temperature value for the vacuum expectation value of the Higgs field 



u = 174 GeV. 
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contributions because of less Boltzmann suppression. On the other hand, the contribution 
from the stop decay is smaller because the process becomes available from the Higgsino- 
axino mixing by the electroweak symmetry breaking and thus involves the top mass rrit 
which is much smaller than fi appearing in the Higgsino decay process. We can observe 
that all contributions become rapidly decreasing after the electroweak scale and drop down 
drastically below 50 GeV being heavily suppressed by the Boltzmann factors. Due to such 
a high Boltzmann suppression, the heavy axino can have a right dark matter density for 
Tr ~ 50 GeV. 



IV. HEAVY AXINO DECAY AND NEUTRALINO DARK MATTER 



In this section, we explore the cosmological impact of a heavy axino which cannot be the 
LSP and thus has to decay to a dark matter candidate which we assume to be a neutralino 
LSP. Such a scenario has been considered first in Ref. 15[ for the case of the KSVZ axino 



and later in Ref. 



12| for the case of the DFSZ axino. Generalizing the latter analysis, we 



examine the dark matter property depending on the neutralino parameter space and the PQ 
scale. 

Depending on the mass spectrum, the axino can decay to a Higgsino and a Higgs, a 
Neutralino (Bino-like) and a pho ton, etc. We assume here, for simplicity and consistency 



with current experiment 



21 



22|, that all the squarks and sleptons are heavier than the 



Higgsino. One can consider two regimes of the axino mass as follows. 



For rria > /i, the dominant decay process of the axino is 



a^H^ + H^/Z or + /W^ 



(8) 



where stands for neutral Higgs states, i.e. /i, if, or A. This decay mode is due 
to the tree-level axino-Higgs-Higgsino coupling. If allowed, it will provide the most 
efficient axino decay. 
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FIG. 8: Contour plot of the axino lifetime (sec). 
For rria < fi, possible decay processes are 

a X + 7 
a ^ Z + X 



(9) 
(10) 

(11) 
(12) 



The first mode is due to the IPI axino-photon-neutralino interaction. The second 
mode is a 3-body decay due to the tree-level axino-Higgs-Higgsino interaction. The 
third mode is a tree-level process due to the axino-Higgs-Higgsino coupling with the 
Higgsino-Bino and axino-Higgsino mixing. Finally the last one is a tree-level process 
due to the axino-Higgsino mixing. 
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Following the calculation described in Appendix, we compute the total decay rate in 
Fig. m Unlike the KSVZ case where the axino yield depends on the reheat termperature, 
the DFSZ axino yield is determined by the mass spectrum and can be parameterized by 

if the reheat temperature is larger than the electroweak scale. Here C, is an order-one quan- 
tity depending on the particle mass spectrum. Such abundant axinos will decay to yield 
overproduction of neutralinos which may thermalize or reannihilate depending on the axino 
lifetime. The axino decay temperature Td is given by 



Td = 1.4 GeV 



If To is larger than the freeze-out temperature Tj of the neutralino LSP, the overproduced 
LSPs will thermalize to be settled down to the usual freeze-out relic density. For To < Tf, 
the overabundant neutra : 
density is determined by 



ino density can be depleted by reannhilation and the final relic 

3 



y^)-^x y^D) + ^^^^^ , (15) 

where Y^{T£)) = y^(Tfr) + ~ 1^ for the reheat temperature above the electroweak scale. 
Since Td < Tf, we need the annihilation rate (cr^f ) of the neutralino LSP larger than the 
canonical value in order to get a right relic density for dark matter. 

For low reheat temperature, Ya can become comparable to or smaller than Y^{T{j.). For 
example, we can find this occurs around the temperature T ~ 40 GeV, 50 GeV, and 60 GeV 
for vpQ = 10^° GeV, 10^^ GeV, and 10^^ GeV, respectively, from the result of Fig. [7] with 
the assumed mass spectrum. Below this temperature, the neutralino relic abudance is given 
by the standard results with a negligible axino contribution. 

In the case of the KSVZ axino, the axino lifetime determined by the coupling ([6]) can be 



15| . On the other 



long enough so that there could occur a period of axino matter domination 
hand, the DFSZ axino lifetime determined by the coupling ([3]) is typically shorter than that 
of the KSVZ axino, and thus there occurs no DFSZ axino domination as shown explicitly in 
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FIG. 9: Black lines stand for the standard freeze-out neutralino abundance. Blue lines stand for 
spin-independent direct detection cross section with proton in unit of pb. The WMAP dark matter 
abundance is depicted by thick line. 

Figs. [T0lfT2l To get an rough estimate, let us consider the largest axino lifetime ~ 10^^ sec 
found in Fig. [I2] for vpq = lO^^ geV. This give Td ~ 0.1 GeV from Eq. dTl). This can be 
compared with the axino-radiation equahty temperature T^g determined by T^g = irriaYa/S 
[isl which give Tgg ~ 10~^ GeV for ~ 1 TeV and Ya ~ 10~^. Thus, it becomes clear to 
have Teg < T^, that is, the axino decays before it overdominates the Universe, for almost 
all the parameter space. Note that the axino lifetime can become much larger near the 
boundary line between the neutralino LSP and axino LSP, so the axino dominated universe 
can be realized. In this case, however, the neutralino abundance becomes too large to get 
right dark matter density. 

Let us now analyze the relic abundance of the neutralino dark matter in the simplified 
MSSM parameter space varying the PQ scale vpg. Only relevant MSSM parameters are 
Ml, M2 and fi in our calculation. 

In Fig. [9l we show the contours of the usual freeze-out relic abundance of the neutralino 



LSP and the direct detection rate in Mi-fi plane, setting M2 = 2 TeV for simp^ 
are numerically obtained by using SUSY-HIT 2J| and micrOMEGAs 2.4 



icity. Contours 



23|. Note that 



the WMAP consistent relic density QuAih'^ ~ 0.1 25[] is obtained in the Higgsino-like LSP 
region (Mi ^ fi) and the resonant annihilation region with Mi ^ mA/2 as expected. One 
can also see that the direct detection experiment [26| excludes the Higgsino-like region and 
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FIG. 10: The left panel shows the neutralino relic abundance after the axino decay, and the right 
panel shows the axino lifetime for vpQ = 10^^ GeV. Red line is the border line between the axino 
LSP and the neutralino LSP. 
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FIG. 11: Same as the previous figure with vpq = 10 GeV 



the lower fi part of the resonance region. 

Let us proceed to consider the effect of the axino late decay in Figs. [TUlfT^ for the reheat 
temperature above the electroweak scale. We take rria = 500 GeV and tan /3 = 10 for all the 
plots, and vpQ = 10^° GeV in Fig. HOl vpQ = 10" GeV in Fig. [H] and vpQ = 10^^ geV in 
Fig.[T2l 

As can be seen in the left panel of Fig. [101 the Higgsino-like LSP with /i ^ 400 GeV 
in the upper- left part gives the correct dark matter relic density. In this part, the Bino- 
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FIG. 12: Same as the previous figure with vpg = 10^^ GeV. 



Higgsino mixing becomes negligible, and the axino decay is determined by the tree-level 
axino-Higgino-Higgs coupling (|3]). This is why the contour is almost parallel to the vertical 
axis in the large Mi region. Note that the phase space of the axino 2-body decay to the 
Higgsino-like LSP and the Higgs becomes very small for fi ^ 400 GeV. For fi > 400 GeV, the 
phase space of the axino 2-body decay is closed, but the axino can decay to the Higgsino-like 
LSP through the processes like 3-body decays with almost massless quarks and leptons via 
the off-shell Higgses or Z-boson, 2-body decay by the axino-photon-neutralino interaction in 
Eqs. ([9]), etc. However all these processes are much suppressed by phase space factors or the 
1-loop factor. So in this region, the axino lives long enough to decay after the Higgsino-like 
LSP freeze-out. Thus, one can understand that the vertical line with fi ^ 400 GeV with 
Ml > 400 GeV is the parameter region which provides a correct dark matter density with 
an appropriately enhanced annihilation rate for the given axino lifetime Tq ~ 10~^ s (as 
shown in the right panel) according to the reannihilation equation ( !T5|) . On the other hand 
one can see that the thick contour line for the correct dark matter density in the lower-right 
part of the left panel is almost same as in Fig. [9] as the corresponding axino lifetime is short 
enough for the axino to decay before the Bino-like LSP freeze-out. 

For vpQ = 10^^ GeV, the axino decays later as seen in the right panel of Fig. [TTl As a 
consequence, the contour line for the right dark matter density changes slightly compared 
to the previous case. More dramatic change can be found for vpQ = 10^^ GeV as in Fig. [T2l 
The Higgsino-like dark matter becomes significantly lighter (/i ~ 250 GeV) and the Bino- 
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like dark matter needs more fine-tuned resonance both of which guarantee more enhanced 
reannihilation by a factor of Tf/Tr,. Note also that the light Higgsino-like dark matter is 
allowed by the direct detection experiment only for the region of Mi > 700 GeV and the 
Bino-like dark matter region requiring stronger resonance is more restricted (compare with 
the blue lines in Fig. 

V. CONCLUSION 

We examined thermal production of the DFSZ axino and its impact on the dark matter 
property. The DFSZ axino interactions are governed by two important couplings: the axino- 
Higgs-Higgsino and axino-top-stop coupling proportional to /i/fpg and mt/vpq, respectively. 
The latter coupling, which arises from the axino-Higgsino mixing, is effective only after the 
electroweak symmetry breaking. All the scattering, decay and inverse decay processes if 
allowed are shown to produce the typical DFSZ axino density of Ya ~ 10~^(10^^GeV/fpQ)^ 
when the reheat temperature is larger than the electroweak scale. In a certain mass param- 
eter region, a resonant scattering can occur to allow enhanced axino production. Thus, if 
the axino is the LSP, its mass has to be very light to get the right dark matter density. If 
the axino LSP mass is at the TeV scale, the reheat temperature needs to be below around 
50 GeV in order not to overdose the Universe. 

A heavy DFSZ axino can decay to a neutralino LSP through the axino-Higgs-Higgsino 
coupling, which leads to a huge neutralino population. Its impact is studied to identify the 
Bino-Higgsino LSP region consistent with the observed dark matter density and the direct 
detection limit of dark matter. For lower PQ scale, a new Higgsino-like LSP region appears 
to meet an appropriate reannnihilation condition, whereas the usual Bino-like LSP region 
with a resonant annihilation through a CP-odd Higgs boson remains unchanged as the axino 
decays before the neutralino freeze-out. For higher PQ scale, the axino decay occurs much 
later than the neutralino freeze-out which requires stronger reannihilation. As a result, the 
neutralino dark matter region is more depleted to a lighter Higgsino-like LSP and a more 
restricted Bino-like LSP. 
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Appendix 



A. Axino-neutralino mixing 



In the supersymmetric DFSZ axion model, we have axino-neutralino mixing because of 
vacuum expectation values of Higgs doublets. Let us write down the axino-neutralino mass 
matrix in the basis of (a, B, W^, H^, H^): 



M 



m-a 





(16) 



( m-„ Cu^VSiijvpQ Cu^VCfilvp(\ 

Ml —CjsSw'mz spsy/mz 

M2 cpcwniz -spcwmz 

CHlJ^vs^/vpQ -cpswmz cpcwmz -fi 

\cH^ivcp/vpQ spswmz -spcwmz -/i / 

where sp = sin/3, cp = cos/3, sw = sin 9w, cw = cos 9w and v = 174 GeV. The above 
matrix can be diagonalized by the method of perturbative diagonalization as in Ref . 19| : 



Mdi^s = NMN'^ = VM'V^ = VUMU^V^ 



(17) 



where 



U 



/l \ 

10 

10 

75 



\0 



(18) 



V2 Vl 



and 



A B 

B^ C 



(19) 
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with 



^ms ^ 



A 



Ml 
M2 



B 



-mzsw{,sp + cp)/\/2 mzsw{sp - cp)/\/2 
mzCw{sp + cp)/V^ mzcw{sfj - cp)/y/2 J 



In the leading order we have 



ML 



-Mnn -MjYim 

As far as \Mi — \ma — /i| ^ v/vpg, fi/vpq, we find 



^3 = 






rria- fi 


V2vpQ{fi - m-a 


Vo4 = 








V2vpq{ij. + rria)' 


Vis = 


M',i 


mzsw{sp + cp) 


Ml- fi 


V2(/i-Mi) ' 


Vu = 


Mil 


mzsw{sp + cp) 


Ml + /i 


V2{fi + Mi) 



where the zeroth component denotes the axino component. Hence we have 



No3 
iVo4 
iVl3 
iVl4 



1 

1 

1 

1 
71 



(K)3 + ^4) = 
hVo3 + ^04) 
{Vi3 + Vu) = 

i-Vis + Vu) 



nit 



mzSwjfJ'Sp + Micp) 

/i2 - Mf 
_ mzSw{MiSp + /icg) 
/i2 - Ml 



C 



fi 

-12^ 

(20) 
(21) 



(22) 
(23) 
(24) 
(25) 



(26) 
(27) 
(28) 
(29) 



Having obtained the mixing matrices N, U and V, one can now get the following axino 
couplings following Ref. [20 1. 

Axino-neutralino-Z boson coupling: 



Cw 



(30) 
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where 

O'L = ^(-A^osiVns + iVo4iV„4). (31) 
Axino-chargino-W^ boson coupling: 

Cw^ax^ = 92W-dY [O^Pl + OIPr] xt + h.c, (32) 



where 



V2 

1 

Axino-neutralino-neutral Higgs coupling: 



O^^ = -^No4Vn2 + No2VnU (33) 



= i=iVo3?7„2 + iVo2f/ni. (34) 



Aax" = 92Thonhdxl, (35) 
^HaxO = g2TH0nHdx^, (36) 
^Adico = mTAOnAd-i^xl., (37) 



where 



T/,o„ = sina(5o„ + cosaSo^, (38) 

Tmn = -cosaQo„ + sina^o^, (39) 

Taou = - sin /3g'o'„ + cos f3Sl, (40) 

Qon = ^iVo3(iVn2 - tanOwN,,,) + i"^^ iVooiVn4 + (0 ^ n), (41) 
^ V2g2VpQ 

S'L = JiVo4(A^n2 - tan ^^^iV^i) - "^"^^ iVooA^ns + (0 O n). (42) 
^ V2g2VpQ 

Axino-chargino-charged Higgs coupling: 

Cn^dx^ = -g2H~h[Q'tPL + Q'oIPr]x^ + h.c, (43) 
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where 



Q;f„ = cos/3[iVo4Ki + -^{No2 + tanewNoi)Vn2] + sin/3^^iVooK2 (44) 

V / 92VPQ 

QZ = sin (3 [NosUni + 4s(iVo2 + tan ewNoi)Un2] + cos /3^^iVoof/„2 (45) 

92VPQ 

B. Axino decay rate 

We first consider the case of rria < /i. Among the 4 possible decay modes, ([9]) is one-loop 
suppressed and ( ITOl) is suppressed by 3 body phase space. On the other hand, (fTTl) and 
( !T2|) are tree-level processes which arise from the axino- Higgsino and Higgsino-Bino mixing. 
Being 2-body decays they are the dominant processes. The decay rates are given by 



m, 



,AV2(i, 4, 4) g2 ^^^^ ^ j(i + _ 411(1 

(46) 



167r k| 



647r 

T^/~ V /, I ~^ ''^aA^/^(l, K^, K^) ,2//. , n2 21 

r(a^/i + x) = ^g2\ThQi\ |(l + /tx) - 

167r 



(47) 



647r vlhr^u? 



T{a^A + x) = ^^^^^^^^^ 



(48) 



64:71 vIqIi"^ 



(49) 



where (72 is SU{2) couphng constant, cw = cosOy/, tw = tan^^y, 9w is the weak mixing 
angle, A^/^(l,a,6) = [(1 - a - b)^ - 4a6]^/^ kz = mz/nia, = M^/rria, = rnh/nia, 
f^H = fTiH/^ai and ka = niA/ma. Ogf , T^oi, Thqi and T401 are given in (155]) . (jSHD and (HUll . 
respectively. The approximation in the second line of each decay rate is valid for <^ 1. 

For the case of > /i, the axino can directly decay to a Higgsino-like neutralino and a 
chargino so that the decay rates can be even higher than those in the case of < The 
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axino decay rates are 



r(5 ^W^ + H^) 



r(a ^H^ + H^) 



m,AV2(l,4,K^) 



gl[\noz? + \THoA']{{l 



IGtt 
327r vIq 

maA^/^(l, Kjj, nf^) \2j_ \rp \21 ( (T j_ ^. \2 ^.2 \ 

1 fl^nia 



(50) 



32vr vIq 



(51) 



maAV^(l,4,^^) 

IGtt 
1 /i^rraa 
327r 



^72'[|^A03p + |T^04p]{(l-«:M)'-4} 



(52) 



2_ ^2\ o 



5 '^Z; ^/iJ 92 



167r 

maAi/2(^l,K^,K2) 

—^2 



[iOo'3p+io^'4n 



„„.{(l + «:,)2-4}{(l-'^.)' + 24} 



(53) 



167r 



{lOo^iP + IOo^iP} 
+ {OoWi}l2«:;. 



{(l-«^^)^ + (l + ^^-2«:^)«:M 



^2 

327r VpQnia 



(54) 



m, 



aAV^(i,4^,4) 2 

1671 ^2 

X {ig^iP + «r}{i 

1 fJ'^rria 



K, 



(55) 



where = mw/ma, i^h± = 'mH±/f^a and = ix/rria- 
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